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Using a chiral unitary approach for meson-baryon scattering in the strangeness zero sector, where 
the N* (1535)Sii resonance is dynamically generated, we study the reactions jp — * n°r)p and jp — > 
■k°K°T, + at photon energies at which the final states are produced close to threshold. Among several 
reaction mechanisms, we find the most important is the excitation of the A* (1700)1)33 state which 
subsequently decays into a pseudoscalar meson and a baryon belonging to the A(1232)P33 decuplet. 
Hence, the reaction provides useful information with which to test current theories of the dynamical 
generation of the low-lying 3/2~ states. The first reaction is shown to lead to sizable cross sections 
, and the iV* (1535)Sii resonance shape is seen clearly in the r/p invariant mass distribution. The 

f^**) ■ same dynamical model is shown to lead to much smaller cross sections at low energies in the second 

reaction. Predictions are made for cross sections and invariant mass distributions which can be 

£S) , compared with ongoing experiments at ELSA. 

C$ • PACS numbers: 25.20.Lj, 11.30.Rd 

(N . 

I. INTRODUCTION 

(N \ 

The unitary extensions of chiral perturbation theory UxPT have brought new light in the study of the meson- 
baryon interaction and have shown that some well known resonances qualify as dynamically generated, or in simpler 
words, they are quasibound states of a meson and a baryon, the properties of which are described in terms of chiral 
■ Lagrangians. After early studies in this direction explaining the A(1405)Soi and the 7V*(1535)Sii as dynamically 
generated resonances 0, 0, 0, 0, , more systematic studies have shown that there are two octets and one singlet 
lO ■ of resonances from the interaction of the octet of pseudoscalar mesons with the octet of stable baryons 0, 0] • The 
A r *(1535)S'ii belongs to one of these two octets and plays an important role in the wN interaction with its coupled 
channels r/N, KA and KT, In spite of the success of the chiral unitary approach in dealing with the meson-baryon 
interaction in these channels, the fact that the quantum numbers of the N*(l535)Sn are compatible with a standard 
three constituent quark structure and that its mass is roughly obtained in many standard quark models 0, , or 
recent lattice gauge calculations nas as a consequence that the case for the N*(1535)Sn to be described as a 
dynamically generated resonance appears less clean than that of the A(1405)iS , u i where both quark models and lattice 
calculations have shown systematic difficulties 01 ■ Ultimately, it will be the ability of the models to describe different 
experiments in which the resonances are produced that will settle the issue of what represents Nature better at a 
certain energy scale. A detailed description of many such experiments has been discussed in ^| . 
' A good example is found in the recent experiment on photoproduction of the A(1405)«S'oi resonance in the 7 p — * 
K + nY, reaction |l4j . where theoretical predictions usin g th e chiral unitary approach had been done previously |l5j|. 
In the present paper we adopt and extend the ideas of and study the analogous reaction -fp — > ir°rip where the 
rjp final state can form the N* (1535)15*11 resonance. This reaction is currently being analyzed at ELSA[la| 

Some of the reaction mechanisms in our model are described as a two-step process: In the initial photoproduction, 
two mesons are generated, one of which is the final 7r°. The final state interaction of the other meson with the proton 
is then responsible for the 77 production. For this interaction chiral Lagrangians in SU(3) representation involving 
only mesons and baryons are used. In addition, the contributions from explicit baryonic resonance exchange such as 
A(1232)P33, A r *(1520)£'i3, and A*n700)Z?33, which have been found essential for the two meson photoproduction, 
e.g., in the Valencia model [H E EH, will be included. The A*(1700)Z? 3 3 resonance, as recent studies show |2fll2lj . 
qualifies as dynamically generated through the interaction of the 0~ meson octet and the 3/2+ baryon decuplet . In 
this picture it is possible [2l| to obtain the coupling of the A*(1700)£> 33 to the 77A(1232)P 33 and if £*(1385)Pi 3 for 
which experimental information does not yet exist. 
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In Sec. [H] the model for the dynamical generation of the AT*(1535)5n resonance in ttN — > ttN scattering will 
be briefly reviewed with special emphasis on the ttN — ► r/N transition. Subsequently, we study the one-meson 
photoproduction -fp — ► r\p in Sec. Ill Al This allows for a simultaneous description of the existing data for these three 
different reactions within the chiral model. In Sec. IHII we predict observables for the photoproduction of ir°r]p in the 
final state. 

At the same time we also study the jp — ► ir°K°'E + reaction and make predictions for its cross section, taking 
advantage of the fact that it appears naturally within the coupled channels formalism of the jp — > ir°r]p reaction and 
leads to a further test of consistency of the ideas explored here. 

In this section we have referred to all resonances that will enter the evaluation of our amplitudes. In what follows 
for shortness of notation we will omit the description in terms of £27,27- 



II. THE iV* (1535) IN MESON-BARYON SCATTERING 

Before turning to the photoproduction reactions of the next sections, let us recall the properties of the A^*(1535) 
in the meson-baryon sector, where this resonance shows up clearly in the spin isospin (S = 1/2,/ = 1/2) channel. 
In the past, this resonance has been proposed to be dynamically generated [ll 0,11,13 rather than being a genuine 
three-quark state. The model of Ref. [g provides an accurate description of the elastic and quasielastic ttN scattering 
in the Sn channel. Within the coupled channel approach in the SU(3) representation of Ref. 0, not only the ttN 
final state is accessible, but also KT,, KA, and rjN in a natural way. 

In the case of the present reactions, we are interested in the rjp interaction which will manifest the A^*(1535) 
resonant character. This interaction was studied in detail in Ref. || for the charge Q = 0, strangeness zero sector. In 
the present study we work in the charge Q = +1 sector, which requires the simultaneous consideration of the coupled 
channels 

TT°p, Tr + n, r}p, K + A, . (1) 

We will subsequently refer to these channels as one through six in the order given above. In this section we derive 
the necessary modifications of the coupled channels in the Q = +1 sector and briefly review the basic formalism. The 
theoretical framework of the photoproduction mechanisms is found in subsequent sections. 

We thus begin with the lowest order chiral Lagrangian for the meson-baryon interaction [22l I23I ] 

C[ B) = {Bi-fVpB) - M B (BB) 

+|d (B 7 "75 K, B}) + \f <B7 m 7s[«m» B \) (2) 

where the symbol () denotes the trace of SU(3) matrices and 

V tl B = d fi B + [T fi ,B], 
1 

2 ( 

U = u 2 = exp(i\/2$//), 

= iu + d^Uu + (3) 



r M = ■^(u + 9 ai u + m5 /j m + ), 



with $ and B the usual 3x3 SU (3) matrices of the fields for the meson octet of the pion and the baryon octet of the 
nucleon, respectively (S^. The term with the covariant derivative V M in Eq. J5J) generates the Weinberg- Tomozawa 
interaction and leads to the lowest order transition amplitude 

Va = -c. l3 -L- u{p')^u{p) (fc M + *;) (4) 

where p,p' (k, k') are the initial and final momenta of the baryons (mesons). The coefficients CV, are SU(3) factors 
which one obtains from the Lagrangian, and the fi are the it, 77, K decay constants |24j . The CV, coefficients for the 
channels with charge +1 are shown in Tab. [I] The amplitudes after unitarization are given in matrix form p| by 
means of the Bethe-Salpeter equation 

T(V5) = [1 - V(Va)G{^] ^ V(V~s) (5) 

with V obtained from Eq. (0J). We are only interested in the s-wave meson-baryon interaction to generate the Sn 
amplitude; the projection of V into this partial wave is given in ref. Q, as well as G which is the meson-baryon loop 
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TABLE I: CVj coefficients for the six channeis. The matrix is symmetric. 
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function in dimensional regularization. In the following we denote by the matrix elements of T with the channel 
ordering of Eq. 

A second modification of the model of Ref. @ with respect to other approaches concerns the ttA^ — > tt-kN channel. 
This channel was important to obtain a good description of the 7 = 3/2 amplitude but it has only a small influence in 
the 1=1/2 channel. It increases the width by about 10 % and changes the position of the A^*(1535) by about 10 MeV. 
In the charge +1 sector this channel can be included by a change of the potential according to V^n^n — * V^m^n + SV 
as in Ref. 123 and reads: 



SV(ir°p -> ir°p) 
6V(n p -> TT+n) 



SV(ir + n — > n+n) 




tA' 



(6) 



with the isospin classification and conventions as in Ref. |8|; G^n being the irirN loop function that incorporates 
the two-pion relative momentum squared. Analytic expressions for Vn and v$i are found in Ref. 

The ttN — > rjN production cross section has been calculated in Ref. and was found to be quantitatively correct 
at the peak position, although somewhat too narrow at higher energies. The question is whether this is due to higher 
partial waves that enter at larger energies and are not part of the calculation, or due to a too narrow iV*(1535) of 
the model. This can now be answered because an 5n partial wave analysis has become available (2(|. In Fig. [TJthis 
analysis is compared to the model of Ref. @- With the solid line, the full model is indicated, and with the dashed 
line the model before introducing the vector exchange in the t-channel and the ttN — * ttttN channel (details in Ref. 
Q)- We will refer to this second one as a "reduced" model of ref. Q in what follows. Although we prefer the full 
model, as form factors and 717^ production certainly play an important role, we take the differences between the 
models in this work as an indication of the theoretical uncertainties. 

In Fig. ^ the energies close to threshold and in particular the strength are well described by the dynamically 
generated resonance. The position of the resonance in the analysis [2(J is at slightly higher energies than predicted 
by the model and the width is considerably larger. This might be due to the contribution of the A r *(1650) resonance 
which is near the iV*(1535) in the Sn channel and has been found to contribute to the reaction in other work [27j . 
Note, however, that in the same reference the total cross section above s 1 / 2 around 1650 MeV is dominated by heavier 
resonances from other partial waves such as the Pi3(1720) and Z?i3(1520). It is also worth noting that in some variants 



4 




1500 1520 1540 1560 1580 1600 1620 1640 1500 1520 1540 1560 1580 1600 1620 1640 



s 1 ' 2 [MeV] s 1 / 2 [MeV] 

FIG. 1: The Sn partial wave in tvN — » r/N. Dots: Analysis from Ref. 12(1 . Solid line: full model from Ref. Dashed line: 
model from Ref. y| without t channel vector exchange and ttttN channel. 




FIG. 2: Photoproduction of rjp via the iV*(1535) resonance (gray blob). Kroll-Ruderman term (a) and the meson pole term 
(b). 

of the chiral models with additional input to the one used here, one can account for the TV* (1650) contribution to the 
Sn amplitude psj. 

In the present approach we restrict ourselves to the model for the TV* (1535) from Ref. |8| as the behavior near the 
rjp threshold is well described in that work including the strength at the maximum of the cross section. 

A. Single meson photoproduction 

In the previous section we have seen that the dynamically generated TV* (1535) resonance provides the correct 
strength in the 7r7V — > rjN transition at low energies. Here, we test the model for the reaction jp — ► rjp with the basic 
photoproduction mechanisms plotted in Fig. [21 which consist of the meson pole term and the Kroll-Ruderman term 
- included for gauge invariance - followed by the rescattcring of the intermediate charged meson described by the 
model of the last section. We shall come back to the question of gauge invariance later on in the section by looking 
at other, subdominant diagrams. 

The baryon-baryon-meson (BBM) vertex is given by the chiral Lagrangian 

£bbm = (Bj^u^B) + (Bj^Bu^) (7) 

with the notation from Sec. [n] The Kroll-Ruderman term is obtained from this interaction by minimal substitution 
and the 7MM couplings emerge from scalar QED. For the ith meson-baryon channel from Eq. the T- matrix 



TABLE II: Isospin coefficients for the Kroll-Ruderman term (a^n,, &kr)i BBM vertex (a BBM , b BBM ), and -yMM vertex (c^ MA/ ). 
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elements read 



V2ie _ / , D + F ul D-F\ , l3) . f d 3 q M 
= j- ^ [aW + V KR J T^(^) J 

?MP{V~S) = ~J~ «( a BBM + ^BBM ^T^) Hw^m)^ (V^) 



(2tt)3 2w(q)£(q) ^fs - E{q) - «(q) + ie 



M t 2 /■ J ^ 2 (l-^ 2 ) 



q /da; 



2(2tt) 2 7 ""7 £(g>) Vi - w - + V« - - k - E{q) + ie 

-1 

^ fc -,-^ + , e ^ + + J) + (w + ^ (8) 

for the Kroll-Ruderman term and the meson pole, respectively. The amplitudes T^ t3 '(y/s ) are the strong transition 
amplitudes from channel i to the rjp channel, following the ordering of table 1 . In Eq. © and throughout this study 
we use the Coulomb gauge (e° = 0, ek = 0, with k the photon three- momentum). The assignment of momenta in Eq. 
© is given according to Fig. [2J u = y/ q 2 + m 2 , a/ = y/ 'q 2 + k 2 — 2qkx + m 2 , E(q) the baryon energy, y/s = P° + fc° 
and G being the meson-baryon loop function according to Ref. The coefficients a, b, c are given in Table [H] The 
cut-off in Eq. JSJ has been chosen A = 1400 MeV. With this value, the reduced model of the rescattering (see comment 
below Eq. (jfJJ)) provides the same strength as the data j2?| at the maximum position of the total cross section. Once 
the cut-off has been fixed, we continue using this value for A in the following sections, for the reduced and full model. 

The amplitudes are unitarized by the coupled channel approach from Ref. [j| in the final state interaction which 
provides at the same time the rj production. This is indicated diagrammatically in Fig. [21 with the gray blob. The 
total amplitude including the rescattering part is then given by 

6 

T 7 „(Vi ) = ) + AtpW*)- (9) 

i=l 

The resulting cross section is plotted in Fig. [21 together with the data compilation from Ref. [^j . With the 
solid line, the result including the full model for the MB — > r/p transition according to Sec. |H]is plotted (dashed 
line: reduced model). The diagrams from Fig. [21 together with the unitarization, explain quantitatively the one- 
meson photoproduction at low energy which indicates that these mechanisms should be included in the two-meson 
photoproduction reactions of the next section. 

Instead of our microscopic description of the r\ production, one can also insert the phenomenological irN — > -qN 
transition amplitude from Sec. HT1 and Ref. (2(| into the rescattering according to Fig. [21 The channels K + YP and 
K + A in the first loop play an important role and should be incorporated as initial states in the MB — > rjp transition. 
In this case we include them by replacing in Eq. © by 

r< a >(V5)-» 5jb^ r (a) (Vi) ( 10 ) 
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FIG. 3: Cross section for 7p — > rjp. Dots: data from Ref. |2flj |. Solid line: Prediction including the full model from Ref. 0. 
Dashed line: Reduced model from Ref. ||. Thin dotted line: Phenomenological ttN — > r)N potential from from Ref. [2(1 . 

where is the phenomenological Sn amplitude to the transition n + n — > r/p. The prescription of Eq. I|1(J|) is the 

correct procedure for the ir + p channel which is the dominant one, and we assume it to be valid for the other channels. 
We choose A = 1400 MeV for the cut-off as before. The cross section is displayed in Fig. [3] with the thin dotted line 
and indeed shows a wider shape. 

As we can see in Fig. [21 the description of the data is only qualitative. Given the theoretical uncertainties one should 
not pretend a better agreement with the data. Yet, in both theoretical calculations the distribution is too narrow, 
reflecting most probably the lack of the N* (1650)Sn contribution in the theoretical calculation. The uncertainties of 
the model for this reaction will be considered later on in the study of the jp — > ir°r]p reaction in order to estimate its 
theoretical uncertainties. 

At this point we would like to make some general comments concerning basic symmetries and the degree to which 
they are respected in our approach, as for instance chiral symmetry or gauge invariance. 

In our approach we are using chiral Lagrangians which are used as the kernel of the Bethe Salpeter equation 
and which arc chiral symmetric up to mass terms which explicitly break the symmetry. The unitarization does not 
break this symmetry of the underlying theory since it is respected in chiral perturbation theory (xPT), and a perfect 
matching with xPT to any order can be obtained with the approach that we use, as shown in Ref. [{j- 

Tests of symmetries can be better done in field theoretical approaches that use, for instance, dimensional regular- 
ization for the loops. Although dimensional regularization is used here in the loops for meson baryon scattering, we 
have preferred to use a cut off for the first loop involving the photon and do some fine tuning to fit the data. Then, 
we use this cut off (which is well within reasonable values) for the other loops that we will find later on. The cut off 
method is also easier and more transparent when dealing with particles with a finite width as it will be our case. The 
use of this cut off scheme or the dimensional regularization are in practice identical, given the matching between the 
two loop functions done in Section 2 of Appendix A of Ref. |30j . There, one finds that the dimensional regularization 
formula and the one with cut off have the same analytical properties (the log-terms) and are numerically equivalent 
for values of the cut off reasonably larger than the on shell momentum of the states of the loop, which is a condition 
respected in our calculations. By fine tuning the subtraction constant in dimensional regularization, or fine tuning 
the cut off, one can make the two expressions identical at one energy and practically equal in a wide range of energies, 
sufficient for studies like the present one. Of particular relevance is the explicit appearance of the log-terms in the 
cut off scheme which preserve all the analytical properties of the scattering amplitude. 

The equivalence of the schemes would also guarantee that gauge invariance is preserved with the cut off scheme if 
it is also the case in dimensional regularization. This of course requires that a full set of Feynman diagrams is chosen 
which guarantees gauge invariance. At this point we can clearly state that the set of diagrams chosen in Fig. [21 is 
not gauge invariant. Some terms are missing, which we describe below, and which are omitted because from previous 
studies we know they are negligible for low energy photons [3ll | . Since the energy of the photon is not so small here, 
it is worth retaking the discussion which we do below. 

The issue of gauge invariance for pairs of interacting particles has received certain attention [32, IU 03 ■ 
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FIG. 4: Photon coupling besides the diagrams from Fig. The photon can also couple to the external baryon (a), internal 
baryon of the first loop (b), and components of the rescattering (c)-(e). 



for the purpose of the present paper we can quote directly the work of |36j| which proves that when using the Bcthe 
Salpctcr equation with the kernel of the Weinberg- Tomozawa term, as we do here, gauge invariance is automatically 
satisfied when the coupling of the photon is made not only to the external legs and vertices, but also to the vertices 
and intermediate particle propagators of the internal structure of the Bethe-Salpeter equation. 

A complete set of diagrams fulfilling gauge invariance requires in addition to the diagrams shown in Fig. [21(a), (b), 
other diagrams where the photon couples to the baryon lines, vertices, or the internal meson lines from rescattering. 
We plot such diagrams in Fig. All of them vanish in the heavy baryon approximation. This is easy to see. In 
diagram (a), Fig. 0] the first loop to the left (think for the moment about a ir + n loop) contains a p-wave vertex of 
the er • q type and an s-wave vertex, and vanishes in any case. Diagram (b) contains a p-wave and an s-wave vertex in 
the loop plus a jnn vertex proportional to a x k. In the baryon propagators one momentum is q and the other one 

q + k and the integral docs not vanish. However, the contribution is of the order (^xr) or ^ nc t crm ( c ) nas 
the same property, a p-wave and an s-wave vertex in the first loop to the left, and only the fact that the propagator 
depends on k + q renders a small contribution (remember we are performing a nonrclativistic calculation by taking 
(7 ■ q for the Yukawa vertices, but this is more than sufficient for the estimates we do). In diagram (d) the MMBBj 
vertex is of the type (<x x q) • e (see Sec. IIII A|l . hence once again we have the same situation as before for the first 
loop to the left. Finally, in diagram (e) the photon is coupled to the internal meson line of the rescattering. In this 
case both the loop of the photon as well as the first one to the left contain just one p-wave coupling and the diagram 
is doubly suppressed. 

In order to know more precisely how small are the diagrams in our particular case we perform the calculation of 
one of them, diagram (b), explicitly. By assuming a 7r + n in the loop to the left in diagram (b), we obtain for the loop 

fi„ r- D + F f d 3 q 



2U J (2tt)3 

1 M M 1 I 

2w(q) E n {q) E N (k + q) - w(q) -k- E N (k + q) + it ^ - w(q) - k - E N (q) + ie Q ' 

while the equivalent loop function for the Kroll-Ruderman term would be 

= -a-eeV2 £±* f *«- - J ^ ^ (12) 

where \i n is the neutron magnetic moment. The explicit evaluation of the terms and f( KR ) indicates that when 
the two terms arc added coherently there is a change of 6% in |i| 2 with respect to the Kroll-Rudcrmann term alone. 
If we add now the term with 7r°p in the intermediate state of the first loop and project over I — 1/2 to match with 
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FIG. 5: Photon interaction with mesons and a baryon. The straight dashed line symbolizes an outgoing meson and the curved 
line the meson in a loop of the final state interaction. 



rjN in the final state, the contribution of the magnetic part is proportional to 2fj, n + fj, p instead of 2^t„ with the 7r + n 
state alone, and the contribution becomes of the order of 2%. The convection term e(p + p') • <t j (2M) (p, p' nucleon 
momenta) of the jpp coupling (not present for the neutron) leads to an equally small contribution. 

The exercise tells us how small is the contribution that vanishes exactly in the heavy baryon limit. Since we do not 
aim at a precision of better than 20 % these terms are negligible for us and hence are not further considered. 



III. ETA PION PHOTOPRODUCTION 



Having reviewed the single rj production in the meson-baryon sector and having applied the model to the single 
r\ photo production we turn now to the more complex reaction jp — > 7r°r)p. The reaction will be discussed in three 
steps: In the first part, the participating hadrons will be only mesons and baryons with their chiral interaction in 
SU(3). In the second part the contributions from explicit baryonic resonances will also be taken into account as they 
are known to play an important role, e.g., in the two pion photoproduction |l7l Il8l ITflj. Finally, the decay channels 
of the A*(1700) into ?7A(1232) and K£*(1385) will be included. 



A. Contact interaction and anomalous magnetic moment 



One of the important features of the models for reactions that produce dynamically generated resonances is that 
the Lagrangians do not involve explicitly the resonance degrees of freedom. Thus, the coupling of photons and mesons 
is due to the more elementary components, in this case the mesons and baryons, which are the building blocks of the 
coupled channels and which lead to the resonance through their interactions. 

We follow the formalism of Ref. 0] for the jp — > K + ttY, reaction where the A(1405) resonance is clearly visible 
in the 7rS invariant mass distribution. The derivative coupling in the meson vertex of Eq. (@J) leads to a ^MMBB 
contact vertex through minimal coupling, sec Fig. [31(c), and guarantees gauge invariance together with the meson 
pole terms of Fig. [3] (a),(b). 

The contact term of Fig. |SJ (c) is easily generated and assuming the reaction jMiBi — ► MjBj the amplitude is 
given by 



,(7) _ r e 

a -°«4A 

with Qi,Qj the meson charges. In the Coulomb gauge this becomes 



Vff' = c„ — (Q, + QMp'h'Mi'K (13) 



y " = - c «iM if ^ f <* +<M (14) 

in the jp CM frame. Since the initial channel i is tt°p, or channel number 1 in the order of the channels from Sec. [HJ 
we obtain 

<"=- a 'S i ¥'' ,i ' (15) 

It was shown in Ref. 0] that the meson pole terms of Fig. [5] (a), (b) are small compared to the amplitude of Eq. 
(|15|l for energies where the final particles are relatively close to threshold, as is the case here, both at the tree level 
or when the photon couples to the mesons within loops. The coupling of the photon to the baryon components was 
also small and will be neglected here, as was done in Ref. |Tfij . 
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Before we proceed to unitarize the amplitude, it is worth looking at the structure of Eq. H15(l which contains the 
ordinary magnetic moment of the proton. It is logical to think that a realistic amplitude should contain also the 
anomalous part of the magnetic moment. This is indeed the case if one considers the effective Lagrangians given in 
Ref. Ill 



4A/ p 



b°(B[S^S»]{F+,B}) (16) 



with 



= d„A v - tf A M (17) 

with M p the proton mass and the electromagnetic field. The operator Q in Eq. (|17fl is the quark charge 
matrix Q = diag(2, —1, — 1)/3 and S 11 is the spin matrix which in the rest frame becomes (0, a/2). In Ref. j3^| the 
Lagrangians of Eq. (|16|) were used to determine the magnetic moment of the A(1405). In the Coulomb gauge one has 
for an incoming photon 

^^^^(Jxq)? (18) 
and, thus, the vertex from the Lagrangian of Eq. (|16f) can be written as 



b% {B{(vtQu + uQiJ) ,#}))• (19) 



+ 2 

Expanding the terms up to two meson fields leads to contact vertices with the same structure as Eq. <|14ll . Taking 
it = 1 in Eq. I|19|) j and hence with no meson fields, provides the full magnetic moments of the octet of baryons from 



where one obtains the values of the coefficients [Hzl I 



= 2.40, b% = 1.82. 

It is easy to see j^l that by setting b$ = 0, b$ = 1, one obtains the ordinary magnetic moments of the baryons 
without the anomalous contribution. Similarly, taking the same values of bg , b§ one obtains Eq. i|14|) for the vertices 
'yMMBB. This is easily seen by explicitly evaluating the matrix elements of Eq. i|19|) which lead to the amplitude 

- ni = - oTtC^ x qK 0x7: + Y ^ ( 20 ) 



where the coefficients Xy and Yij are given in Table UTTI The combination of the Y\j in Table UTTl and the C\j of Table 
HI shows the identity of Eq. J5D|) and Eq. (pjj for the case of b® = 0, bg =1. 



1. Unitarization 

For the amplitude 7p — > ir°iip, the first thing to realize is that at tree level the amplitude is zero with the interactions 
from Eqs. (|15fl an d H2(J|) . It is the unitarization and the coupled channel procedure that renders this amplitude finite 
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FIG. 6: Unitarization of the transition amplitude for r)p production. The possible states (see table^J for one and two loops are 
indicated. 



and sizable. The unitarization procedure with the coupled channels allows the intermediate channels with charged 
mesons of Eq. Q to be formed, even if some of them are not physically open. The scattering of these states leads 
finally to i]p. Diagrammatically, this is depicted in Fig. which implicitly assumes the unitarization is implemented 
via the use of the Bethe-Salpeter equation JSJ) which generates the diagrams of Fig. 
Since r\p is channel 3 in our list of coupled channels, our final amplitude reads 

iwv = + r«) j^t ' Gj(z) tU3){z) ' (21) 

where Gj is the meson-baryon loop function which is obtained in Ref. @ using dimensional regularization, and the 
are the ordinary scattering matrices of the rjp and coupled channels from Eq. ©. The invariant kinematical 
argument z is given by the invariant mass Mi{rjp) of the r\p system, 

z = Mi (22) 

or, alternatively, 

z=(s + ml-2V- S pl) 1/2 (23) 

with p° = (p^ + m 2 ) 1 / 2 , when the amplitude is expressed in terms of the invariant mass Mi(n°p) of the n°p system. 

One might also question why we do not unitarize the other ir° with the r\ or the proton. The reason has to do with 
the chosen kinematics. By being close to threshold the ir° has a small momentum and is far from the region of the 
ao(980) resonance that could be created interacting with the rj. The generation of the ir°p invariant masses in the 
A(1232) region in the phase space that we investigate is more likely. However, as one can see from Fig. 5 an extra 
loop of the 7T° and p lines produces a A(1232) which would involve an s-wave vertex and a p-wave vertex. This would 
vanish in the loop integration in the limit of large baryon masses. Later, we shall consider other diagrams in which 
the A(1232) is explicitly produced. 



B. Kroll-Ruderman and meson pole term 

Next, we take into account diagrams which involve the 7^V — > tjN amplitude which has been discussed in Sec. Ill Al 
and which are shown in Fig. [7\ With Eqs. © and J^J the amplitude for the diagrams in Fig. [3 is given by 

D + F M i 
( ^ = ^(k + Pw ) E N (k)-pO-E N (k + P7T ) (-* T -rP^vp(z)) ' P.) (24) 

where z takes the values given by Eq. 12 I'll or (|23|l . 



1. Intermediate pion emission 

In addition to the diagrams considered above, there are additional diagrams in which the 7r° is produced inside the 
first meson-baryon loop as displayed in Fig. [H] The amplitude for the channel i for the sum of diagram (c) and (d) is 
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FIG. 7: Kroll-Ruderman term and meson pole term as sub-processes in the n°ri production. 





'7T°, Ptt '7T , Ptt 



I I 

FIG. 8: Pion emission from inside the first meson-baryon loop. Diagram (d) is required by gauge invariance. 



given by 



t 



(c) + (d) 



(a ■ P7r )(ae) [a t (D + F) + b t (D - F)][a{ (D + F) + b[ (D - F)] 



d 3 q 1 



M 



M' 



(2tt) 3 2lo E M (q) 
1 



■p° - £ M '(q + pJ + «e 



1 - 



E M (q) ~ 

gon 



ie£? M '(q + P w ) 
F n (q-k) 



(25) 



where the index i stands for our standard ordering of the channels in Eq. and the only non-zero values of the 
a l ,a! i ,b i ,b' i are: a 2 = ~^i a 2 = - l> a 4 = b 4 = ^,04 = v/| , b' A = --^,a 5 = 65 = ^,65 = ^775- N °t e tnat channel 

two has the external 7r° coupled to n, n to the left and right in the diagram, channel four has the 7r° coupled to the 
A, S to the left and right, and channel five has the 7r° coupled to the E, A to the left and right. In the equation the 
variable P — q refers to the baryon on the left (M) of the emitted ir° and the variable P — q—p n to the right (M 1 ) of 
the emitted tt° as shown in Fig. [S] The contribution of the terms in Fig. |HJ is therefore given by the sum of Eq. I|25|) 
for the three non-vanishing channels. 

In Eq. 112511 we introduce the ordinary meson-baryon form factor F^ of monopole type with A= 1.25 GeV as used 
in the two pion photoproduction |l7j . It appears naturally in the meson pole term of Fig. [S]and, as done in Ref. 
|l7j . it is also included in Fig. [Sfc) (Kroll-Ruderman term) for reasons of gauge invariance. This form factor does not 
change much the results and it is approximated by taking the q° variable on shell and making an angle average of the 
q momentum. This is done to avoid fictitious poles in the q° integrations. 

The meson pole term (d) in Fig. [S] is small and an approximation can be made for the intermediate pion at 
q — k which is far off-shell. This concerns terms with mixed scalar products of the form k • q that give only a small 
contribution when integrating over q in Eq. (|25|l . Additionally, we have set in this term q = q on , the on-shell 
momentum of the other meson at q. This is, considering the kinematics, a good approximation. 

One can also have the pion emission from the final proton. However, this would imply having the 7rA^ — > rjN 
amplitude away from the A^*(1535) resonance at a value Mi = ^fs where the 7rAf — > rjN amplitude would only 
provide a background term above the A^*(1535) resonance. Once again the set of diagrams considered leads to 
small cross sections compared to the dominant terms to be considered later in the paper so further refinements arc 
unnecessary. 
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FIG. 9: Terms with A*(1700), JV*(1520), and A(1232). The diagram on the right is the A-Kroll-Ruderman term. The latter 
diagram implies also a meson pole contribution, required by gauge invariance, which is not separately plotted but is included 
in the calculation. 



C. Baryonic resonances in 7771-° production 

In the present study, the 7771" production is described as a two-step process: The first step consists in the photo- 
production of two mesons and a baryon; the second step describes the subsequent transitions of meson-baryon — > ijp 
via the dynamically generated iV*(1535) resonance. In particular, the first stage contains two pion photoproduction. 
For this part it is known that baryonic resonances such as A's and N*'s can play an important role [T7L ITsL ITfll. Isfij. 
For this reason wc include the relevant mechanisms from Ref. [T^J adapted to the present context. Fig. [5] shows the 
processes that are taken into account. The s-wave character of the AT* (1535) (gray blob in Fig. [§)) discards all those 
processes from Ref. |T^ where both pions couple in p-wavc to the baryons, because the irN loop function involving 
odd powers of q in the integral is zero in the heavy baryon limit. For the remaining processes, some contributions to 
the 7T + 7r° and 7r°7r° cross sections are small as, e.g., from the Roper resonance. Finally, one is left with the A*(1700)A, 
AT* (1520) A, and A-Kroll-Ruderman terms from Fig. [UJ The latter implies also a pole term which is required by gauge 
invariance in the same way as in Fig. |Hl(d). Since many resonances appear in this section we refer the reader to the 
notation used in the Introduction. 

The amplitudes for diagrams (e) and (f) in Fig. [5] are given by 

A 1 

i — 1,2 n 1 



2 



with the sum only over the first two channels according to Eq. Q and A = 1400 MeV as in Sec. Ill Al The meson 
energy ui, baryon energy E and energy of the A(1232) read 

2 2,2 
to = m + q , 

E 2 = M 2 +cf+pl +2qp^x, 

sa = (v 7 ^-^) 2 -Q 2 (27) 

with meson mass m, baryon mass M, and p w = \p T \ with p\ = \Jp\ + m 2 . The argument z is given by Eq. (|22(l or 
(Ho- 
using the notation from Ref. [l7| the amplitudes t l A , which can depend on the loop momentum, are given by: 

£ A-(1700) +t N*(l520) +t A-KR (28) 
,2 _ t A*(1700) t JV*(1520) , ^A— KR ^9) 
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where 



A*(1700) 
~7P — >7r+7r°n 



t 



2 /. 



ANir 



<h 



2M 



1 ffA'A^ 

3 ml 



(axk)e-i&-e(g' 1 (k° + ^) 



N* (1520) 
'7p— >7r+7r°n 



t 



. V2 f. 



= -i ^ S -pn[f 



3 m 



N*'At\ 



1 9n'' At\ 
3 mi 



f 



A-KR 
7P — ^7r+7r°n 

A* (1700) 
7P — >7r 7r n 

TV* (1520) 



51 2M 

9 
1 



(axk)e-i&-e( gi (k + —) 



AN-n 



l 



(30) 



(2p,r - i(a x p w )) • e F^g^ - fc) 1 - 



9on 



jp — *7r°7r°n 



2V2 

V2t 



^A*(1700) 

7P — >7r + "7r°n ' 



iV*(1520) 

+ 7r°n' 



7P — >7r~ r 7r u n' 



A-KR 



= 0. 



(31) 

(32) 

(33) 

(34) 
(35) 

We have already projected out the s-wave parts of the A*(1700)A7r and 7V*(1520)A7r transitions that come from the 
term .9/v*'Air/ m iT ^ ' V-k S ■ p w , see Ref. |l7j . The vector p^ depends implicitly on the invariant mass which will be 
specified later, Eqs. I|46[) or l|48|) . The amplitudes in Eqs. (|30[1 - H35|) are formulated for real photons, which is the case 
we are considering here. The meson pole diagram related to the A-Kroll-Ruderman term has been included in the 
last factor of Eq. (|32|l by making the same approximation as in Eq. (|25|l for the intermediate off-shell pion. The pion 
form factor (see Ref. 01) nas to be inserted since the intermediate pion in the meson pole term is far off-shell. 
For the A* propagator, 

Ga- Wi ) = - - 1 r(V - s y (36) 



Ma* 



Ait (Nitty) we obtain in a similar way as in Ref. |l7j 

qCM{y/s ) 5 



the (momentum-dependent) width according to its main decay channels has been taken into account: For A* — > TtN 
in d-wavc, A* — > Np(N-rnr) 1 and A* 

Fa^NttIVs ) = rA»->iV7r(MA*) 
^A*^N p [tttt](Vs ) 

A = 



M N m A * 
6(2rr) 3 ^ 



g 2 f 2 / dm dui 2 \D p ( qi + q 2 )\\c ll - q 2 ) 2 6(l - \A\), 



(■sfs-wx- uj 2 ) 2 - M 2 N - q? - q| 



rA*^A7r[Ar7T7r] 



15 

16^ 



dMi 



2|qi||q 2 | 

Mi fc(Mj) r A ^AT^(A//) (\A S \ 2 



{Mi - Ma) 2 + ( rA - w 2 ' (Mj) ) 



0(yfs — Mi - m-x) . 



(37) 

Here, 9cm(v / *) * s the CM momentum of the pion and the nucleon and T a*^n-k{M&*) is determined through the 
branching ratio into that channel. For the decay into Np, g p = 2.6 is the A*Np coupling, also determined through 
the branching ratio. Furthermore, /„ = 6.14 is the p7T7r coupling, qi = (w^q^), i = 1,2 the four-momentum of the 
outgoing pions, and D p the p propagator incorporating the p width. For the decay into A7r, the finite width of the A, 
rA-+JV7r, has been taken into account by perf orming the convolution. For the partial amplitudes A s and Ad of the A* 
decay into A7r in s and <i-wave, see Ref. |l7|. The iV*(1520) propagator is dressed in a similar way with the analytic 
expressions given in Ref. [l7j . 



D. 



5(7(3) couplings of the A* (1700) and background terms 



In Ref. |2l| the rescattering of the meson octet with the 3/2+ baryon decuplet leads to a set of dynamically 
generated resonances, one of which has been identified with the A* (1700). The advantage of such a microscopic model 
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FIG. 10: Left side: Coupling of the dynamically generated A* (1700) to KY,* and r)A. The loop is given by rjA + p, K + T,*°A, 
or _?«: E* + E + . Right side: E* Kroll-Ruderman term. 



is that couplings of the resonance to decay channels are predicted which have not yet been determined experimentally. 
In particular, the analytic continuation of the amplitude to the complex plane provides at the pole position the isospin 
3/2 couplings of the resonance to r]A and isTE*. Identifying the pole with the A*(1700) we can incorporate the model 
from Ref. [2l[ in the present study in the diagrammatic way as indicated on the left side of Fig. with the jpA* 
coupling from Ref. |l7j . This procedure can be regarded as a first step towards the incorporation of dynamically 
generated 3/2~ resonances in the two meson photoproduction. In further studies, the initial jp — > A* process could 
be included in the microscopic model of Ref. |2lj in a similar way as was done here for the 7piV*(1535) coupling in 
Sec. Ill Al However, phenomcnologically, the procedure followed here is reliable. 

The A*(1700)Ar/ and A*(1700)T,*K couplings from Ref. H2 are given up to a global sign by g n = 1.7 - HA and 
gx = 3.3 + i0.7, respectively. However, in Ref. [2lJ the coupling to Att is also given, g& = 0.5 + i 0.8. The sign of 
the real part and the order of magnitude agree with the empirical analysis of the A* (1700) — > irA decay that we are 
using thus far |l7j : hence, we take for g n and gx the values quoted above. We note that the cross section is almost 
independent of the global sign, whereas there are some minor differences in the invariant mass spectra. 

Having included the S* in the A* decay it is straightforward to consider also the corresponding S*-Kroll-Ruderman 
term given on the right side of Fig. ^| This term, together with the other ones from this section, allows for an 
extension of the model to higher energies, where the intermediate A(1232) from the processes of Sec. IIII CI is off-shell 
but the £*(1385) is on-shell. 

For the baryon decuplet baryon octet meson octet vertices, and the corresponding Kroll-Ruderman vertex, we take 
the effective Lagrangian from Ref. |4jj j 



C = C 



(38) 



with the same phase phase conventions for the states of the decuplet as taken there, which is the same one taken in 
Ref. [2i|. This allows us to relate all the couplings to the one of AnN. Up to a different phase, these factors agree 
with those used in Ref. |4lj . 

The corresponding amplitudes for the diagrams in Fig. 1101 read now: 



(3) 



(5) 
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(6) 



f (5) 

l E*-KR 



lg n lhllL GA , {V - s) S. Pn 



2D + F 
5 2U 



gi< Ga'(Vs)S ■ p w 



5* -k 



24 D + F 
L15 V ox 9k -777 — G A * {Vs)S ■ p T 
V 25 2/ w 



2M 



(a x k)e - 5t • e(g[(k° + ^j) + 9^ k°^j 



2M (-k).f-5t.d#» + -) + s ^^ 



4V3fD + F\ 2 , ^ ^ 

1-15 e -^r- I — — — I (2p T - 1(0 x p w )) • e 



25 



2/ w 



(39) 
,(40) 
(41) 
(42) 



with g v and gx given in Ref. |2ll |. In order to obtain the full amplitudes, T~ p _n r o r , p , these t's have to be inserted as 
t\ Eq. I|26(l but the sum over index i goes now from three to six. The lower index for the amplitudes in Eqs. (|39() - 
(I42|l indicates the particles in the loop to be considered in the evaluation of Eq. I|26|) . The upper index indicates the 
channel number i and therefore which T^ t3 > has to be chosen in Eq. (|26fl . For the amplitudes from Eq. I|4U[I. Q41JI. 
and l]42p. the A(1232) propagator in Eq. I|2ti|) has to be replaced with the S*(1385) one. The latter is defined in 
the same way as the A(1232) propagator and we take a momentum-dependent width with r rest = 36 MeV assuming 
the dominant p-wave decay of the S* into pA. The numerical factor of 1.15 appearing in Eqs. (|40|) and l|42ll is a 
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FIG. 11: Tree level process from the decay of the A*(1700) to 77A(1232). 



phcnomcnological correction factor from the SU(3) £*7rA coupling in order to provide the empirical E* — > irA partial 
decay width. 

For the ir°r]p production, the A* At? coupling together with the subsequent A — > ir°p decay provides also a term at 
tree level as shown in Fig. 1111 The contribution for this reaction is simply given by 

Tj£L„o np = tfl +p Ga(z') (43) 
from Eq. 1)39(1 . The invariant argument z' for this amplitude amplitude differs from z of the former processes, 



z = [s + m 



,-2^) , z'=Mj (44) 



with p® = (pj? + m 2 .) 1 ' 2 depending on whether the amplitude is parametrized in terms of Mj(r]p) or Mi(n°p), 

respectively. We have explicitly tested that recoil corrections for the A(1232) decay S ■ p v in Eq. I|43|) . in the way 
they are applied in Ref. [13, are negligible. 



IV. RESULTS 



In this section, invariant mass spectra Mj(rjp) and Mj(tt°p) for the reaction ^p — > ir r]p are predicted, together 
with the total cross section for this reaction. The corresponding observables for the tt°K°'S + final state are also given. 
These observables can be directly compared to ongoing experiments at the ELS A facility [lfi|. 

We evaluate the phase space integrals for the invariant mass distribution of r\p in the rjp CM system, 



(1(7 



1 



M„Mi 



2jv 



2tt 



P^Ptt 



dMi(r)p) 4(2tt) 



Ml 



dcj) w / dcosl 



dcost 



(45) 



with p v the modulus of the momentum p^ of the r\ in the r\p rest frame p^ — ^^(Mj ,171^, M^)/ (2Mj) in terms of 

the ordinary Kallen function where the direction of p v is given by <p and 8. This vector is connected to p v in the jp 
rest frame by the boost 



Pr, 



Mr 



Mr 



(-Ptt) +Pr, 



(46) 



where 0^ 



and the tt° three momentum in the jp CM frame is given by the modulus p^ 



\ 1 / 2 {s,Mf,m^ T )/{2yfs) and the two angles <t> n ,0 n . Furthermore, lo^ is the pion energy in the jp CM frame. In 
Eq. I)45|l . M p ,Mi are proton mass and mass of the final baryon, in the present case also a proton (i = 3 with the 
channel ordering from Eq. Eqs. (|45|l and (|46|l arc a generalization of the corresponding expression in Ref. 

as in the present case the amplitude depends explicitly on the angles of the particles. 

The individual numerical contributions from the various processes from Sec. Mil are shown in Figs. 1121 1131 and 
ITll We have chosen here a lab energy for the photon of E~ ( =1.2 GeV so that the allowed invariant mass range is 
wide enough to distinguish the A r *(1535) from pure phase space. On the other hand, this energy is low enough, so 
that unknown contributions from heavier resonances than the A* (1700) should be small. All contributions contain 
the resonant structure of the Af*(1535) in the final state interaction, except the background term from Eq. I(43|l . 
Although the shape of this contribution is similar to the resonant part, this is a combined effect of phase space and 
the intermediate A(1232) that becomes less off-shell at lower invariant masses for E 1 = 1.2 GeV. 
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FIG. 12: Invariant mass at E 1 — 1.2 GeV. Dotted line: Contact interaction from Fig. HJ including the anomalous magnetic 
moment. Dashed dotted line: Meson pole plus Kroll-Ruderman term from Fig. [7] Double dashed dotted line: A*(1700)A'E* 
transitions from Fig. 1101 see Eqs. (1401 . 141H . Solid line: Intermediate pion emission from Fig. |S| 
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FIG. 13: Invariant mass at E-, — 1.2 GeV. Processes with explicit resonances. Dotted line: A*(1700)-7rA contribution from Fig. 
GJ(e), see Eqs. £H3 and J33"). Solid line: Contribution from A^*(1520)?rA in Fig. 0(e), see Eqs. ^ and GHJl. Dashed dotted 
line: A-Kroll-Ruderman term from Fig. Q3(f), see Eq. 1321 . Double dashed dotted line: E*-Kroll-Ruderman term from Fig. 
mil see Eq. lEHZl . 



The individual contributions shown in Figs. El El and 1141 are evaluated using the full model for the iV*(1535) 
from Sec. |H]with the coherent sum indicated in Fig. El solid line. The coherent sum using the reduced model is 
displayed with the dashed line in Fig. El We take the difference between the two curves as an indication of the 
theoretical uncertainty as in the previous sections. 

The first thing to note is that the peak position of the V*(1535) is lowered by some 20 MeV due to the interference 
of the dynamically generated resonance with the background term from Fig. El A width of 93 MeV for the iV*(1535) 
has been extracted in Ref. In the invariant mass spectra the iV*(1535) exhibits a considerably smaller width. 
This is for two reasons: First, the JV*(1535) is situated close to the rjp threshold and the phase space cuts the lower 
energy tail. This is clearly visible in Fig. El As the phase space factors in Eq. I|45|l are smooth functions around the 
iV*(1535) resonance, the shape of the curves in Figs. El El and 1141 reflects the 7V*(1535) resonance seen through a 
\T\ 2 matrix involving the coupled channels. 

The second reason for the narrow iV*(1535) is that at higher invariant mass the amplitude for the resonance is 
suppressed by the initial photoproduction mechanism: A closer inspection of the dominant resonant contributions as, 
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FIG. 14: Invariant mass at _E 7 = 1.2 GeV. Dashed dotted line: A*(1700)r;A transition and r\p — > rjp rescattering from Fig. 1101 
sec Eq. dSHJ. Dotted line: Tree level process with A*(1700)?7A transition but no rescattering, Fig. 1111 see Eq. I|43(l . Solid line: 
Coherent sum of all contributions, full model for the iV*(1535). Dashed line: Coherent sum of all contributions, reduced model 
for iY*(1535) from Sec. |H](no vector particles in t-channel, no -k-kN channel). 




FIG. 15: Modulus of the amplitude (full model for iV*(1535)) for the reaction MiBi -> rjp (left side) and MiBi -t ir°p (right 
side). As these amplitudes appear squared in invariant mass spectra and total cross sections, they serve as a useful tool to 
distinguish dominant processes for the ty tjp final state. Initial states: Dotted lines: Dashed lines: iv + n, Solid lines: rjp, 

Dashed dotted line: K + T,°, Double dashed dotted lines: K + A, Triple dashed dotted lines: K°T, + . 



e.g., from Eq. I|39[l shows that the A(1232) propagator in Eq. (|2(j[l of the first loop becomes more and more off-shell 
at higher rjp invariant masses which leads to a suppression of the spectrum for this kinematics. This effect is in fact so 
pronounced that the shape of the invariant mass distribution hardly changes if the MB — > r\p theoretical transition 
amplitudes are replaced in the scattering diagrams by the phcnomcnological ones given by Eq. (|l(Jf) . This is clearly 
seen in Fig. ^| 

The transitions |T( 83 )| and \T^ \ in Fig. [^explain the size of some of the contributions in Figs. and 1 141 as 

they appear squared in invariant mass spectra and cross section. E.g., the ttN — > rjp transitions on the left side of Fig. 
E3are small which explains why the A-Kroll-Ruderman term and the 7V*(1525)7rA, A*(1700)7tA transitions from 
Eqs. H3l)fl - 135f) contribute little, opposite to what was found in the two-pion photoproduction [17j|. In contrast, the 
diagrams using SU(3) Lagrangians without explicit resonances from Figs. EH3 and |5] contain KA and KY, channels 
in the first loop so that the contributions are larger. The by far largest contributions in the rescattering part (Fig. 
1141 comes from the A* (1700) — > 77A decay with the subsequent unitarization of r\p. Indeed, the r/p — > rjp scattering 
amplitude is very large as Fig. 1151 shows. Additionally, the loop for this reaction in Fig. 1101 contains a 77 instead of 
a 7r, and the particles in the loop can be simultaneously on-shcll, whereas for the 7rAA^ loop at least one particle is 
always further off-shell. 

The diagrams with E*(1385) in the first loop are relatively large (Fig. EJ> due to the large A*(1700)A'E* coupling 
and the large KY, — ► rjp and KA — ► r\p transitions from Fig. ED However, the S*(1385) is off-shell at E 1 = 1.2 GeV 
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FIG. 16: Phenomenological potential for the MB — > r\p transition at E-, = 1.2 GeV. With solid and dashed line, full and 
reduced model for the iV*(1535) as in Fig. 1141 Dotted line: Phenomenological potential in the meson-baryon — > ryp final state 
interaction for the diagrams from Figs. EJUHHHID and llOl 




FIG. 17: Selection of diagrams with n p being the final state of the rescattering instead of rjp. These diagrams are suppressed. 

and the contribution can not become as big as the loop from the A* (1700) — > Arj decay. Therefore, diagrams with 
a A*(1700)A"£* coupling become more important at higher energies. From Fig. right side, we can also directly 
read off that additional diagrams like those displayed in Fig. 1171 which use instead of are small compared 
to their counterparts from Sec. IIIII 

A. Extension to higher energies 

In Fig. I18l the results for the invariant mass distribution arc shown for higher values of the incoming 7 momentum, 
?iab = 1-2 — 1.7 GeV. The resonant shape of the TV* (1535) is not modified if a bigger photon energy is chosen, only the 
size decreases slightly as the intermediate A* (1700) of the dominant processes becomes off-shell. At higher incident 
photon energies, the peak of the iV*(1535) moves back to its original position around 1520-1540 MeV (see, e.g., Fig. 
^) as the interference of the dynamically generated iV*(1535) with the tree level process from Fig. Illlbecomes weaker. 

A second maximum appears for E~ t >~ 1.5 GeV and moves to higher invariant masses with increasing photon energy. 
This can be traced back to be a reflection of the A(1232) resonance in the tree level process from Fig. El which is 
on-shell around the position of the second peak. When predicting this double hump structure in the r/p invariant 
mass, one has to keep in mind that our model for the dynamically generated A r *(1535) resonance underpredicts the 
width of this resonance (see, e.g., Figs. ^13- Furthermore, there are unknown contributions from resonances heavier 
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FIG. 18: Invariant mass spectrum 
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and dashed lines: Full and reduced model for the TV* (1535), respectively. 
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FIG. 19: Integrated cross section for the 7p — » 7r°?7p reaction. Solid line: Full model for the V*(1535). Dashed line: Reduced 
model (see Sec. |HJ. Dashed dotted line: Phenomenological potential for the MB — > rjp transition (only available up to E~, ~ 1.2 
GeV). 



than the A*(1700) about which little is known and which can fill up the space in invariant mass between the two 
humps. As a result, we expect a separation of the two maxima not at E 7 = 1.5 GeV as Fig. 1181 suggests but at higher 
energies. Nevertheless, the tree level process from Fig. EJcontributes so strongly to the coherent sum that the double 
hump structure should be qualitatively visible in experiment. 

In Fig. ^5]the integrated cross section is shown. There is a steep rise below E 1 — 1.2 GeV simply due to growing 
phase space. Above that, the cross section grows slower and finally saturates. At high photon energies, the tree level 
process and the dynamically generated iV*(1535) are almost completely separated in invariant mass (see Fig. I18f) 
and we do not expect a further rise beyond 1.7 GeV within our approach, as the particles involved in the various 
processes become more and more off-shell. However, the narrow A^*(1535) width of our model, together with unknown 
contributions from resonances heavier than the A* (1700), lead to uncertainties at high photon energies which are hard 
to control. 

As we have already seen in Fig. ^j] the use of the wider phenomenological potential increases the cross sec- 
tion slightly (dashed dotted line), but it is remarkable how insensitive the cross section is to the actual width of 
the iV*(1535), regarding the large difference in width between the results using the phenomenological potential or 
microscopic theory which we have seen for the 7p — > rjp reaction in Fig. [3J 
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FIG. 20: Invariant mass spectrum dMj d ^o p ) [m d GeV 1 ] as a function of M/(7r°p) [MeV] for various photon lab energies _E 7 . 
Solid lines: Full model for the 7V*(1535). Gray lines: Phase space only (T=const). Dashed dotted lines: Only tree level process 
from Fig. 1111 see Eq. (03- Dotted line in plot for E^ — 1.2 GeV: Effect when including recoil corrections for the A7r°p vertex 
for the tree level diagram in Fig. 1111 



B. The tv°p invariant mass 

In the discussion of the last section we have seen that the A(1232) plays a prominent role in the TT°r]p photopro- 
duction. For completeness, the invariant mass spectra for the 7r°p particle pair is given, which should show a signal 
of the A(1232). The phase space integrals are evaluated in the ir°p rest frame and lead — similar to the expression 
in Eq. gSl — to the invariant mass distribution for Mi(tt p): 



da 



1 



M p M, 



PnPr, 



dM^p) 4(2tt) £ 



Ml 



d cos 9 r 



dcos8j2J2\ T - 



7P — >rjTT u p I 



(47) 



with Mi = M p and with p n the modulus of the momentum 
X 1 / 2 (Mf,ml,Mp)/(2M I ) where the direction of p„ is given by (f> 
7P rest frame by the boost 



p v of the 7T° in the ir°p rest frame p n = 
and 9. This vector is connected to pV in the 



- 1 



PirPv 



Mi 



(-Pv)+P7r 



(48) 



where p° 



P-K 



mi 



and the rj three momentum in the 7P CM frame is given by the modulus p v 



iT)/(2s/s) and the two angles </>,,, 6* l; . Note that the invariant arguments for the solution of the Bcthe- 



Salpeter equation (j3J) have changed compared to the case when the amplitude is expressed in terms of the rjp invariant 
mass, see Eqs. and 

The invariant mass distribution including all processes from this study is plotted in Fig. [201 We have checked 
explicitly for the individual processes and for the coherent sum of all processes that the integration over Mj(tt°p) in 
Eq. I|47|) leads to the same values for the cross section as when integrating over the r\p invariant mass distribution 
from Eq. (|45|l . In the plot for E 1 = 1.2 GeV, the dotted line indicates the negligible effect of recoil corrections for 
the tree level process from Fig. ^2 as described below Eq. (|44|) . 

In Fig. [2n]we observe at E 1 = 1.2 GeV a shift of strength towards higher invariant masses compared to the pure 
phase space (gray line) obtained by setting T = const in Eq. I|47|) . This is caused by the low energy tail of the 
A(1232) from the tree level process from Fig. ^2 indicated with the dashed-dotted line. Indeed, at higher photon 
energies, the intermediate A(1232) in this process shows up as a shoulder at E 1 = 1.5 GeV, and as a clear peak 
beyond. Additionally, there is a shift of strength towards higher invariant masses that results in a maximum which 
moves with energy, as it becomes apparent at E 1 — 1.5 GeV. This is a reflection of the iV*(1535) resonance that 
becomes on-shell around these invariant masses, in full analogy to the reflection of the A(1232) resonance in the r\p 



invariant mass spectra in Fig. 1181 As we have already argued in Sec. IIV Al the separation of the two peaks might 
happen at higher values of the incident photon energy but should be qualitatively visible in experiment. 

At this point we would like to make some comments concerning the accuracy of our results. If one looks at the 
results obtained for the jp — > rjp cross section in Fig. [21 we can see that except in the low energy regime close to 
threshold, wc have large discrepancies between the three options and also with experiment. The agreement can be 
considered just qualitatively. It is clear that some background and the contribution of the iV*(1650)5n might be 
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FIG. 21: Tree level contribution for the ty K°T, + final state. 



£ 7 = 1.5 GeV 





1700 1720 1740 1760 1780 

M/(i4: E+) [MeV] 



1700 1725 1750 1775 1800 1825 1850 1875 

i\/ / (A' E+) [MeV] 



FIG. 22: Invariant mass spectrum for the reaction jp — > 7r A"°E + for two photon lab energies E- 1 . Solid lines: full model for 
the iV*(1535). Dashed lines: reduced model for the A r *(1535). Dotted lines: Only tree level process from Fig. 



missing and that in any case the theoretical uncertainties from different acceptable options are as big as 20-25 % at 
some energies. We should not expect better agreement with experiment in the jp — ► n ijp reaction which requires the 
IP ~ * VP amplitude in some terms (see Fig. J7J. However, as we have discussed, these terms give a small contribution 
to the total amplitude, since the largest contribution comes from the tree level diagram of Fig. Illl and its unitarization 
in Fig. 1101 Thus, at the end, the uncertainties in the result for the r\p invariant mass distribution, as seen in Fig. 
1161 are smaller than those of Fig. |3| Furthermore, when one integrates over the r\p invariant mass distribution, the 
uncertainties in the calculation in the total cross section are relatively small, although we would not claim a precision 
of better than 20 % considering all the different sources that enter the calculation. Given the complexity of the model, 
such an uncertainty is not easy to decrease at the present time, but it is more than acceptable for this first model of 
the reaction. 



C. The reaction jp -t tt°K ^ 



The 7p — + tt°K°'E + reaction is calculated in a similar way as in the last sections for the Tr°r]p final state as the 
coupled channel formalism for the iV*(1535) contains the K°Y, + final state in a natural way. There is, however, a 
different tree level diagram as displayed in Fig. [21] with the amplitude 
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(49) 



(50) 



in analogy to Eq. 144|) . when expressing the amplitude in terms of the K°Y, + invariant mass. The £*(1385) propagator 



Gs* has been given its width as explained below Eq. (|42|) 



Note that T BG „ „ + 



from Eq. (|41|l in analogy to Eq. (|43|l . For the contributions with rescattering, we simply have to choose the (i,6) 
channel instead of the (i,3) channel in from Eq. (JSJ, in the ordering of the channels from Eq. This means 
the replacement T^ 3 ' — > T^ 6 * 1 in Eq. I|21|l and accordingly for the rest of the contributions. The invariant mass 
distribution is obtained from a similar formula as Eq. I|45|) with Mj = Af s + and is plotted in Fig. 1221 

The tree level contribution from Fig. 1211 dotted line, dominates the spectrum. The reaction is situated at much 
higher energies than the jp — > ir°r]p process and the dynamically generated iV*(1535) is off-shell for the entire invariant 



22 



0.3 
0.25 
„ 0.2 

^= 0.15 

b 

0.1 
0.05 


1350 1400 1450 1500 1550 1600 1650 1700 

E 1 [MeV] 

FIG. 23: Integrated cross section for the 7p — > tv°K T, + reaction. Solid line: Full model for the JV*(1535). Dashed line: 
Reduced model for the JV*(1535). Dotted line: Contribution from the tree level diagram in Fig. 1211 



mass range. Thus, the rescattering part appears as a uniform background. The tree level process shows a pronounced 
maximum which moves with the incident photon energy. This situation is analogous to the moving peak of the 
A(1232) in Fig. 1 181 and reflects the £*(1385) which is on-shcll around the peak position. The full and reduced model 
for the final state interaction (see Sec. HI) differ considerably in Fig. El at E 1 = 1.7 GcV (solid versus dashed line). 
This is due to the fact that the model for the iV*(1535) becomes uncertain at these high energies as it differs from the 
ttN — > nN partial wave analysis from Ref. [2i| above y/s ~ 1600 McV. At lower energies, the differences are smaller. 
In any case, in the energy range studied here, the dominant term is provided by the tree level diagram of Fig. 1211 

The integrated cross section is displayed in Fig. [23] Comparing with Fig. ^]it becomes obvious that the production 
of tt°K°Y 1 + is highly suppressed. This is a combined effect of the A* (1700) and the iV*(1535) being off shell which we 
quantify below. Both reactions are compared at an energy of 150 MeV above their respective thresholds, where both 
cross sections have become significantly different from zero. For the resulting photon energies of E 1 = 1202 MeV and 
1603 McV for the ir°rip and tt°K°^ + final states, respectively, we obtain c 7r o ?w ,/er 7r o i< -o s + = 10.9. First, this ratio is an 
effect of the positive interference between the dynamically generated iV*(1535) with the large contribution of the tree 
level term from Fig. ^2 Calculating the cross section by using this latter term only, (7^0,^(1202 McV) decreases by 
a factor 2.0. Second, and more important, the A* (1700) propagator from Eq. I)36|) is off shell for the higher photon 
energy, |G A *(1202 MeV) | 2 /|G A * (1603 MeV)| 2 = 5.4. Multiplying these two factors, one obtains 10.8 which clarifies 
the origin of the factor 10.9 quoted above. 

Turning the argument around, if the experiment sees a factor 10 suppression of the ir°K°'E + final state, compared 
to 7r°7/p, this can be easily explained by the dominant role of the A* (1700) found in the present study. 

Other resonances beyond the A* (1700) can contribute at these high energies, and their omission produces uncer- 
tainties in the calculated cross section. However, assuming that we have included the relevant mechanisms in the 
present model, the suppression of the ir°K°Y, + versus the ir°i]p final state is such a strong effect that it should be 
visible in experiment. 

V. CONCLUSIONS 

In this paper we have studied the reactions jp — > ir°r]p and jp — > ir°K°'E + making use of a chiral unitary framework 
which considers the interaction of mesons and baryons in coupled channels and dynamically generates the TV* (1535). 
This resonance appears from the s-wave rescattering of r\ N and coupled channels. We have used general chiral 
Lagrangians for the photoproduction mechanisms and have shown that even if at tree level the amplitudes for these 
reactions are zero, the unitarization in coupled channels renders the cross sections finite by coupling the photon to 
intermediate charged meson channels that lead to the r\p and in the final state through multiple scattering of 

the coupled channels. 

The theoretical framework has been complemented by other ingredients, considering explicit excitation of reso- 
nances, whose couplings to photons are taken from experiment. 

The interaction of the meson octet with the baryon decuplet leads to a set of dynamically generated resonances, 
one of which has been identified with the A* (1700). The decay of this resonance into rjA and KY,* , followed by the 
unitarization, or in other words, the A* (1700) — » ir°N* (1535) decay, provides in fact the dominant contribution to 
the iV*(1535) peak in the invariant mass spectrum. A similar term provides also a tree level process which leads, 
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together with the iV*(1535), to a characteristic double hump structure in the r/p and ir°p invariant mass at higher 
photon energies. 

A virtue of this approach, concerning the r/p spectrum around the iV*(1535), and a test of the nature of this 
resonance as a dynamically generated object, is that one can make predictions about cross sections for the production 
of the resonance without introducing the resonance explicitly into the formalism; only its components in the (0~, l/2 + ) 
and (0~, 3/2 + ) meson-baryon base are what matters, together with the coupling of the photons to these components 
and their interaction in a coupled channel formalism. The reactions studied here also probe decay channels of the 
A*(1700) -> ??A(1232), A*(1700) -> if E* (1385) or transitions like r/p -> 7V*(1535) ->■ rjp which are predicted by the 
model and not measured yet. 

We have made predictions for the cross sections and for invariant mass distributions in the case of the 7p — » 7r°?7p 
reaction. For the second reaction under study, 7p — > ir°K°Y, + , we could see that in the regions not too far from 
threshold of the 7p — > ir°K°Y> + reaction, the cross section for the latter one was much smaller than for the first 
reaction. 

The measurement of both cross sections is being performed at the ELSA/Bonn Laboratory and hence the predictions 
are both interesting and opportune and can help us gain a better insight in the nature of some resonances, particularly 
the AT*(1535) and the A* (1700) in the present case. 
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